Abstract. We have carried out a quantitative trend anal- 
ysis of the crystahine siHcates observed in the ISO spectra 
of a sample of 14 stars with different evolutionary back- 
grounds. We have modeled the spectra using a simple dust 
radiative transfer model and have correlated the results 
with other known parameters. We confirm the abundance 
difference of the crystalline silicates in disk and in outflow 
sources, as found by Molster et al. (1999a). We found some 
evidence that the enstatite over forsterite abundance ra- 
tio differs, it is slightly higher in the outflow sources with 
. respect to the disk sources. It is clear that more data is 
^? ' required to fully test this hypothesis. We show that the 
. 69.0 micron feature, attributed to forsterite, may be a 
' very suitable temperature indicator. We found that the 
^ . enstatite is more abundant than forsterite in almost all 
l«— ^ ' sources. The temperature of the enstatite grains is about 
. equal to that of the forsterite grains in the disk sources but 
T— I ' slightly lower in the outflow sources. Crystalline silicates 
, are on average colder than amorphous silicates. This may 
' be due to the difference in Fe content of both materials. 
Finally we find an indication that the ratio of ortho to 
clino enstatite, which is about 1:1 in disk sources, shifts 
towards ortho enstatite in the high luminosity (outflow) 
T-H ' sources. 
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1. Introduction 

The Infrared Space Observatory (ISO) has provided a new 
and unprecedented view on the occurence and compo- 
sition of circumstellar and interstellar dust. One of the 
surprises of the ISO mission was the discovery of ubiq- 
uitous crystalline silicates in circumstellar dust shells of 
both evolved and young stars (see e.g. Waters et al. 1996; 
Waelkens et al. 1996). We have carried out an extensive 
study of the presence and properties of crystalline sili- 
cates. The present study is the third in a series, in which 
we study these silicates using ISO spectra. In previous pa- 
pers (Molster et al. 2001c; 2001d; hereafter Papers I and 
II respectively) we have measured and described the cir- 
cumstellar dust features found in the infrared spectra of 17 
stars with different evolutionary status. The majority of 
these features could be identified with crystalline olivines 
(Mg2xFe2-2xSi04) and pyroxenes (Mg:cFei_xSi03), where 
1 > cc > 0. Jager et al. (1998, hereafter JMD) measured 
the mass absorption coefficient of crystalline pyroxenes 
and olivines with different Fe over Mg ratios. Bands of 
both materials show a shift in the wavelength position of 
the peaks to longer wavelengths with increasing Fe con- 
tent. The detection of the 69 micron feature, which is very 
sensitive to the Fe/Mg ratio (Koike et al. 1993; JMD), as 
well as the relative strength of the crystalline silicate fea- 
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tures in the spectrum of IRAS09425-6040 (Molster et al. 
1999a), led to the conclusion that the crystalline olivines 
observed in the ISO spectra are very Mg-rich {x > 0.95); 
the Mg-rich end member of the olivines is called forsterite. 
Similarly, the enstatite band at 40.5/im is sensitive to the 
Fe/(Fe-f Mg) ratio and points top the presence of Mg-rich 
pyroxenes. The identification of the dust species is very 
important for a better insight in the formation and evolu- 
tion of dust. This may lead to a better understanding of 
the mass loss process and thus the evolution of the mass- 
losing star itself. 

There is a clear separation between sources with and 
without a dusty disk. This difference is evident quanti- 
tatively in the sense that the crystalline silicate features 
are stronger with respect to the continuum in the sources 
which are surrounded by a disk (Molster et al. 1999a), and 
also qualitatively in the shape of the features, which is a 
proof for different dust properties (Paper I). 

However, more quantitative statements are necessary 
to come to a better understanding of the nature of the 
circumstellar dust in these objects, and of their forma- 
tion and processing history. In order to get these quanti- 
tative statements, a comparison with laboratory measure- 
ments is necessary. Unfortunately the laboratory measure- 
ments do not always agree with each other. In Paper II 
we discussed different laboratory measurements of olivines 
and pyroxenes, and possible causes of discrepancy. Despite 
these differences, qualitative agreement with the ISO spec- 
tra is already quite impressive as we will demonstrate in 
the present study. 

In Section^ we discuss trends in the position and width 
of the solid state bands. In Section |^ we apply a very sim- 
ple optically thin dust model to the spectra to determine 
a typical temperature for the dust species. The results of 
this modelling are used to look for correlations which are 
discussed in Section 0. 
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Fig. 1. The peak position of the forsterite peak at 23.7 
versus the peak position of the forsterite peak at 33.6 pirn 
(left side) and the forsterite peak position at 69.0 ^im ver- 
sus the position at 33.6 fj,m (right side). The diamonds denote 
the disk sources and the circles denote the outflow sources. 
The open triangles in the left part are room temperature lab- 
oratory observations for crystalline olivines with and 6% of 
[FeO], Mg2Si04 and Mgi.88Feo.i2Si04 respectively. In the right 
part the laboratory measurements fall off scale due to temper- 
ature effects (see Fig. ^ for these measurements). The error 
bars denote 1 a errors in the wavelength position. No obvious 
correlations are visible. 



Peak positions show variation from source to source. 
We adopt the naming in Paper I, which implies that if we 
refer to a wavelength position we use fim, while if we refer 
to the name of specific feature we will write 'micron'. 

2. Peak positions and bandwidths 

As the measurements in paper I show, there is a spread 
in the peak positions and bandwidths of the different fea- 
tures (see e.g. Fig. ^ where we plotted the spread in the 
2 strongest olivine peaks at 23.6 and 33.6 /im). There is 
no clear trend in the spread. We note that, apart from 
89 Her, all disk sources show the 33.6 micron features at 
a longer wavelength than the outflow sources (see Fig. |l|). 

Here we will discuss possible causes for this spread. 
One of the best investigated causes for shifts in bands is 
the chemical composition. Difference in the mineralogical 
composition as well as the abundance of the different el- 
ements are known to change the peak positions. Another 
cause for a change in the peak position is the tempera- 
ture: a lower temperature tends to narrow the features, 
shift them to shorter wavelengths and in some cases even 
split the band. A third method why bands shift and change 
width is variations in the size and shape of the dust grains. 
Changing the crystallinity of the material can also change 
the appearance of the feature. Below we will discuss these 
four effects more extensively. 



2.1. Composition and temperature 

JMD and Koike et al. (1993) showed that the inclusion of 
Fe in (crystalline) silicates will increase the wavelengths 
of the peak positions of the different features. Since the 
shift in wavenumbers is rather constant for the different 
features and proportional to the [FeO] content (JMD), 
these shifts are best seen for the features at the longest 
wavelengths (see Fig. ^). However, even a plot of the 69.0 
micron feature versus the 33.6 micron feature does not 
show a clear trend (see Fig. I) . Both Fig. I and Fig. | show 
that the crystalline olivines are very Mg-rich and Fe-poor. 
The wavelength positions are consistent with an absence 
of Fe in these crystals. We conclude that the spread in the 
observations cannot (only) be explained by a very small 
and changing [FeO] content, therefore another mechanism 
should also be responsible. 

Another mechanism to shift bands is the temperature. 
The 69.0 micron forsterite feature is found in different lab- 
oratory spectra where it always peaks at 69.7 /im, while 
in our ISO spectra it is always found around 69.0 ^m. 
This shift is significant. Lowering the temperature will 
shift this feature bluewards (e.g. Bowey et al. 2000, Chi- 
hara et al. 2001). Another effect of a temperature decrease 
is a narrowing of this feature (see Fig. ||). This is a gen- 
eral property of crystalline features reflecting the anhar- 
monic interaction of the phonon modes with the thermal 
phonon bath. At higher temperature the phonon modes 
are more excited and their distribution is broader. Hence, 
phonon-assisted absorption will shift bands redwards and 
will broaden their profiles at higher temperature. The 
amount of broadening depends on the origin of the feature. 
The 69.0 micron band is one of the best isolated bands in 
our spectrum to test for this effect. This narrowing of the 
absorption bands with decreasing temperature may be re- 
sponsible for the fact that the observed band widths (of 
dust with typical temperatures of « lOOK) in almost all 
cases is smaller than the laboratory widths measured at 
room temperature. 

Besides shifting and narrowing a band, bands can also 
split in two components upon lowering the temperature 
(e.g. Bowey et al. 2000, Chihara et al. 2001). This might 
explain why in some sources we see a blend and in oth- 
ers we find two separate components. The narrowing (and 
splitting) of the features might provide an independent 
measurement of the temperature of the dust, without 
knowing anything of the rest of the spectrum! 

We conclude from these comparisons that the crys- 
talline silicates are very Mg-rich and cold. However, both 
processes shift the peaks together along the same line in 
the wavelength versus FWHM diagram, so other processes 
must play a role too. 
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Fig. 2. The observed FWHM and peak wavelength of the 69.0 micron feature in the spectra of the dust around stars 
(open diamonds for the sources with a disk and open circles for the sources without a disk) and in the laboratory 
at different temperatures (filled triangles - forsterite (Foioo; Bowey et al., 2000), and filled squares - olivine (Fogo; 
Mennella et al. 1998)). The temperatures are indicated at each point, and within the resolution the 24 K and 100 K 
for Fogo are similar. Note that the measurements were not corrected for the instrumental FWHM (w 0.29, for the ISO 
observations, and 0.25 and l.O/im for the laboratory observations of respectively Foioo and Fogo) 



2.2. Shape and size 

Strong transitions as found in the crystalline silicates can 
be very shape dependent. Not only the band strength is 
affected by grain shape, also the peak position can change 
dramatically. For instance, the peak wavelength of the 
33.6 micron forsterite feature is located around 32.7 /j,m 
for spherical grains, while it shifts to 33.8 /zm for a conti- 
nous distribution of ellipsoids (see Fig. H). The same figure 
also clearly demonstrates that the relative strength of the 
features is very shape dependent. 

Another shape effect, which can play a role in the 
shape of the spectrum, is the preferential growth in the 
direction of 1 or 2 crystallographic axes (see e.g. Bradley 
et al. 1983 and reference in there). Although this is not 
likely to shift and broaden features that much, unless fea- 
tures corresponding to 2 different axes blend significantly, 
it does play a role in the relative strength of the features. 

Difference in the size of the particles can also influence 
the spectra. However, this only becomes important when 



the size of the particles is comparable to the wavelength 
of the feature. Relatively large particles will broaden the 
feature and shift it to longer wavelengths. 

Finally, we would like to note that coagulation and 
porousity can also have an effect on the width of the 
bands. This has been theoretically investigated for amor- 
phous quartz spheres by Bohren and Huffman (1983) and 
in the laboratory by Koike and Shibai (1994). This effect 
is difficult to quantify, but is likely to play a role for the 
crystalline silicates. 

To conclude, we can say that shape effects might play 
a role in the spread of the peak positions. 



2.3. Crystallinity 

A final consideration should be the crystallinity, i.e. de- 
gree of lattice order or number density of lattice defects. 
Single crystals without defects show much sharper peaks 
than those with some defects. Therefore the width of the 
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Fig. 3. The absorptivity (Qabs) of forsterite measured in 
the laboratory and calculated from the optical constants 
of Servoin & Piriou (1973) for spherical particles with a 
radius of 0.01 /xm, using Mie calculations, and for a con- 
tinuous distribution of ellipsoids (having a volume equal 
to that of a spherical particle with a radius of 0.01 /im). 
All curves are normalized to the peak value of the 33.6 
micron feature. Note the difference in peak position and 
strength between the 2 different shape distributions. 

peaks might be an indication of the crystallinity of the 
dust grains. 

Fig. ^ shows a comparison between forsterite in differ- 
ent environments and 2 laboratory spectra. The forsterite 
measured by JMD was formed from a melt and proba- 
bly polycrystalline, while the forsterite measured by Koike 
et al. (2000b) was a single crystal. The poly-crystalline 
forsterite is expected to have defects where the different 
crystals meet each other. 

The comparison with the astrophysical spectra sug- 
gests that around evolved stars nature produces nice single 
crystals, while around young stars poly-crystalline mate- 
rial seems to be formed. Around young stars it is expected 
that crystalline silicates form at the inner edge of accretion 
disks where the dust particles will (partially) melt. When 
mixed with cooler regions the liquid drops will likely solid- 
ify again and crystallize. In the outflows of evolved stars, 
gas will slowly cool and condense into dust grains. The 
time and temperature in these environments are likely to 
be sufficient to completely crystallize and get rid of all 
defects. It seems that this is not the case for young stars, 
where the (partially) melted dust grains cool more rapidly, 
being unable to remove all the defects. 

Fig. ^ already shows that the crystallinity might shift 
the peak position of the feature a few tenth of a micron. 
So it is likely to play a role in the spread. 

2.4- Summary 

Above we discussed several mechanisms to shift and 
broaden features. We can conclude that the crystalline 
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Fig. 4. The 33 micron complex of the young star 
HD100546 (Malfait et al. 1998), the comet Hale Bopp 
(Crovisier et al. 1997) and an average of 33 micron com- 
plexes of evolved stars with evidence for a disk (Molster et 
al. submitted to A&A), together with 2 laboratory mea- 
surements of forsterite one done by JMD and one by Koike 
et al. (2000b). Note the difference in width of the 33.6 mi- 
cron feature. 



silicates are Fe-poor and cold. However, the reason for 
the spread in peak position is not well known. Shape ef- 
fects and crystallinity, due to the different environments 
in which these particles are formed, are likely to play a 
role, but a combination with temperature and composi- 
tion effects cannot be excluded. 

Finally, it should be noted that all the features, of 
which the peak positions are plotted in Fig. || are part 
of a complex. Contributions of nearby features from other 
materials, which were not detected as separate features, 
might result in (apparent) differences of the peak posi- 
tions. In this respect it is interesting to note that the 
33.6 micron feature is the dominant feature in the 33 mi- 
cron complex and also seems to show the least spread in 
the observations. The 23.6 and 69.0 micron features are 
part of the 23 and 60 micron complexes, respectively, and 
are much less dominant in these complexes. So relatively 
minor contributions will have less impact on the 33.6 mi- 
cron feature than on the 23.6 and 69.0 micron features. 
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2.5. Other claimed trends 

In a very limited sample Voors (1999) found a constant 
separation (in fim) between the 30.6 and the 32.8 micron 
band, which suggest a common species for these two fea- 
tures. We could not confirm this, but found a relation 
with Xx indicating the wavelength of the x micron fea- 
ture. Although we found a trend for those two features, 
we consider it unlikely that the 2 dust features come from 
the same material, since their strength normalized to the 
33.6 micron feature does not correlate. Although we did 
not further investigate the correlation found, it may reflect 
a common cause, e.g. a difference in general temperature 
of the dust particles. 

Based on a study of the crystalline silicate features 
in AFGL4106, Molster et al. (1999b) concluded that the 
forsterite features were broader than the enstatite fea- 
tures. We have tested whether this conclusion holds for 
our larger sample, using the mean results from paper II, 
and could not confirm their claim. A similar negative re- 
sult was found, when we compared the different laboratory 
datasets. 

3. The modelling 

In this section we will apply a simple dust emission model 
to derive the typical temperature and abundance ratio 
of forsterite and enstatite in the dust shells. We assume 
that the dust shell is optically thin at infrared wavelengths 
(which is reasonable because we see the dust features in 
emission) , the grain size distribution of the different dust 
species is similar to what has been measured in the labo- 
ratory (^ 2/im, which also is quite realistic, because the 
width of the features indicates grains smaller than the 
wavelength) and that all grains of a given composition 
have the same (single) temperature (this is probably less 
realistic, but we only want to get a typical temperature 
for the dust species and are at the moment not interested 
in the temperature distribution). A comparison of sev- 
eral laboratory data sets with the observations indicates 
that the laboratory data of Koike et al. (2000b) give a 
good qualitative match to the observations. We will use 
this data set for our modelling. We determined an eye-ball 
spline-fit continuum, maximizing the continuum and still 
be smooth (no sudden changes in the slope), both in F,^ 
and F\. This continuum was derived in a similar way (us- 
ing similar wavelength positions as continuum) for both 
the stellar and the laboratory spectra. Whenever possi- 
ble we tried to use the whole wavelength range (SWS -I- 
LWS) to determine the placement of the continuum for 
the stellar spectra (see also Paper I). We have fitted the 
continuum subtracted spectra with the continuum sub- 
tracted forsterite and enstatite (50% ortho-enstatite and 
50% clino-enstatite) mass absorption coefficients multi- 
plied by blackbody functions. 

^(i^)inodol = 



-F'(i^)modci is the calculated model flux, B{Ti,v) is the 
blackbody function at temperature Ti of dust species i, 
K,(i')i is the mass absorption coefficient of dust species i 
and Mi is a multiplication factor which is related to the 
total mass of dust species i. 

The temperature of the blackbody is not necessarily 
the same for enstatite and forsterite. In determining the 
best fit, we varied the temperature in steps of 5 K. The 
resulting spectra were separately scaled to fit the spec- 
trum. This scaling factor is related to the mass of the 
dust species. The absolute masses requires knowledge of 
the distances to the stars but, for each source the masses of 
the different dust components can be directly compared. 
The mineral mass ratios determined in this paper assume 
that they have the same grain size and shape distribution 
(both around stars and in the laboratory samples). The 
best fits were determined by eye and no method has 
been applied. This method is of sufficient accuracy given 
the current quality of the lab data and given the fact that 
several prominent dust features still lack identification, 
thus strongly affecting any method. We found that the 
temperature and mass for forsterite could be determined 
using the 23 and 33 micron complexes, while the enstatite 
values are mainly based on the 28 and 40 micron features. 

The results of this simple fitting procedure are shown 
in Fig. ^ to |l6| and the derived temperatures and mass ra- 
tios are given in Table We also derived an estimate for 
the typical temperature of the underlying continuum. For 
this we assumed that the continuum is caused by small 
grains with optical constants based on the amorphous sil- 
icate set 1 of Ossenkopf et al. (1992) and a continuous dis- 
tribution of ellipsoids as shape distribution. We fitted the 
continuum to the original, not the continuum subtracted, 
spectra. An independent fit based on a Q{\) ~ emis- 
sivity law gave similar temperatures. This gave us confi- 
dence that the continuum temperature is reasonably well 
determined in this way. It should be noted that other 
shape distributions (e.g. spheres) and other sets of optical 
constants of amorphous olivines can easily change the de- 
rived temperature by ± 20 K, more often to higher than 
to lower temperatures. From these fits we could in princi- 
ple derive a relative mass, like in the case of enstatite and 
forsterite. Although the uncertainties in the (mass) ab- 
sorption coefficients (due to shape, size and compositional 
differences) are systematic, the spread in values makes it 
very difficult to interpret them and to compare them with 
other observations. Therefore, we have not given an amor- 
phous over crystalline silicate mass ratio. However, since 
the differences between the different datasets are system- 
atic, trends can still be derived from these numbers. For 
the remainder of this paper we will take the tempera- 
ture derived by the fit with the Ossenkopf data set as 
the continuum temperature (Table |I|), because these fits 
tend to produce the best fits. We compared the temper- 
atures found in this study with those found by Molster 
et al. (1999b; 2001a), and we found a reasonable agree- 
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Table 1. The derived temperatures for forsterite Tf, enstatite Tc and the amorphous silicates Ta from our model fits. Also the 
forsterite to enstatite ratio is given. The '-' denotes that it was not possible to derive a realistic value, f indicates that it is not 
really possible to fit the spectrum with a single temperature, The temperatures here are found by a fit to the short wavelength 
side of the spectrum. The typical error in the temperature of the crystalline silicates is 10 K, for the continuum temperature 
about 20 K, and in the mass ratio a factor 2. (P)PN — (proto-)planetary nebula, RSG = red supergiant. 



Star 


Type 






Ta 


Me 
Mf 


Disk sources 


IRAS09425-6040 


C-star with O-dust 


85 


100 


145 


1.2 


NGC6537 


hot PN 


75 


65 


80t 


5.8 


NGC6302 


hot PN 


65 


70 


80^ 


1.0 


MWC922 


PPN?, Herbig star? 


90 


100 


140 


2.7 


AC Her 


RV Tauri star 


100 


90 


225''' 


5.0 


HD45677 


Herbig star? 


140 


140 


235^ 


5.3 


89 Her 


PPN 


110 


100 


320+ 


5.0 


MWC300 


PPN? 


90 


90? 


145 


1.4? 


HD44179 


C-PPN with O-dust 


135 


135 


120 


4.0 


Outflow sources 


HD161796 


PPN 


105 


80 


100 


11.4 


HD179821 


post-RSG? 


75 


65 


90 


3.3 


AFGL4106 


post-RSG 


100 


80 


120 


8.0 


NML Cyg 


RSG 


150 




180 




IRC+10420 


post-RSG 


90 




160 





ment. Difference in the temperatures found could often 
be described to the use of different laboratory data sets. 

Our simple model, consisting of only two crystalline 
dust components and a single temperature for each dust 
component, fits most stars very well, see e.g. MWC922 
(Fig. H). Still, it is clear that this simple model is not 
sufficient to explain all the features. The main discrep- 
ancies between our model fits and the ISO data lie at 
the wavelengths below 20 iim. We note that the three 
stars with a continuum temperature above 200 K all show 
crystalline silicates in emission in the 10 micron region. 
The temperature of the crystalline silicates has been de- 
termined based on bands at wavelengths longwards of 20 
micron. These bands are dominated by cool dust, and the 
derived low temperatures (Table ^ are too low to explain 
the strength of the crystalline silicate bands in the 10 mi- 
cron complex. A second, much warmer, component must 
be introduced to explain these 10 micron bands. Likely 
a temperature gradient is present in these sources. The 
discrepancies shortwards of 20/xm do not solely reflect the 
presence of a temperature gradient in these sources, but 
indicate that still other dust components must be present. 
The 18 micron complex is badly fitted. The modelled 19.5 
micron feature (originating from both forsterite and en- 
statite) is often much too strong and the modelled 18.0 
and 18.9 micron features are often too weak when com- 
pared to the ISO spectra. The too strong 19.5 micron fea- 
ture might be a radiative transfer effect, since this feature 
is less of a problem in the full radiative transfer modelling 
(see e.g. Molster et al. 1999b; 2001a). This might indicate 
that our assumption of r ^ 1 is not correct at wavelengths 



around 19.5 /im. The poor fit of the 18.0 and 18.9 micron 
features suggests the presence of another dust component. 

There is more evidence for the presence of an extra 
dust component. The 29.6 and 30.6 micron features also 
need extra emissivity, as is very clear in the spectra of 
NGC6537 (Fig. ^ and of NGC6302 (Fig. 0). In these two 
sources the 40.5 micron feature is not well fitted, suggest- 
ing that the same dust component which is responsible for 
the 29.6 and 30.6 micron features also has a peak around 
40.5 ^m. A possible candidate for this extra dust com- 
ponent is diopside (MgCaSi206), which peaks at the re- 
quired wavelengths. However this material also produces 
strong peaks at other wavelengths, e.g. at 20.6, 25.1 and 
32.1 /im, which are observed in the ISO data, but often 
not as strong as expected. Therefore, the identification of 
the carrier of the 29.6 and 30.6 micron features remains 
open. It should be noted, that the temperature and rela- 
tive mass of enstatite are estimated from the 28 and 40 
micron complexes. Therefore a significant contribution of 
an unknown dust component to one (or both) of these 
2 complexes can change the estimated temperature and 
abundance of enstatite. 

The 33.0 micron feature is not well fitted, but this 
feature is likely to be influenced by instrumental behaviour 
(see Paper II). In the 35 micron plateau we clearly miss 
intensity around 34.8 /im in all sources. The predicted 69.0 
micron feature is often too weak with respect to the ISO 
spectra (see e.g. Fig. |5|). This may be an indication for 
the presence of colder dust, and thus for a temperature 
gradient, we will come back to this later. 

Apart from all the features that are missing, we also 
have a problem with too much intensity predicted by our 
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Fig. 5. A fit (dotted line) to the continuum subtracted spec- pig. 6. A fit (dotted line) to the continuum subtracted spec- 
trum (solid line) of IRAS09425-6040. = 85 K and = trum (solid line) of NGC6537. T, = 75 K and = 65 K. 
100 K. 



modelling around 27 /im. This excess is mainly due to 
enstatite, but also forsterite contributes slightly. We are 
still looking for an explanation of this phenomenon. 

Finally, we did not attempt to fit the absorption pro- 
files. As stated above we assumed the dust was optically 
thin. Also, no attempt was done to fit the carbon dust 
features, which are present in some sources. 

3.1. The sample stars 

Here we describe the model fits to the spectra of the in- 
dividual stars, which where analyzed in Paper I. For a 
description of the individual stars in this sample we refer 
to Paper I. From this sample we rejected Roberts 22 and 
VY 2-2, because the ISO satellite was unfortunately offset 
when observing these two objects resulting in large flux 
jumps in the overall spectrum. This made it impossible 
to derive temperature estimates of the dust around these 
two stars. OH26.5+0.6 has also not been fitted, because 
below 30 /xm it has an absorption spectrum (Sylvester et 
al. 1999), which could not be described with our simple 
model. 

The main uncertainties in the model fits are due to 
uncertainties in the continuum subtraction. This leads to 
errors in the temperature of the order of 10 K and mass 
uncertainties of the order of a factor 2. We note that for 
our modelling we completely rely on the laboratory data 
input. This may result in systematic effects on our derived 
temperatures and masses. 

3.1.1. IRAS09425-6050 

The fit to the spectrum of IRAS09425-6040 is shown in 
Fig. Ij. The model fit also produces a somewhat too strong 
19.5 micron feature. It should be noted that the full ra- 
diative transfer calculations of Molster et al. (2001a) pro- 



duces excellent fits to the 19.5 micron feature. The broad 
feature at 11 /im is due to SiC. This very simple model 
predicts no significant flux in the 10 micron complex due 
to crystalline silicates, which is consistent with its absence 
in the ISO spectrum. 

The forsterite grains have a temperature of 85 K. This 
temperature agrees with the temperature range presented 
in the detailed radiative transfer model (Molster et al. 
2001a). However, in contrast to the results presented here 
these detailed calculations predict that enstatite is much 
cooler than forsterite. As a result those models could not 
reproduce the relative strength of the observed 28 and 
40 micron complexes. It also resulted in an unrealistically 
high mass for the enstatite. Molster et al (2001a) argue 
that this might have to do with the not well known ab- 
sorptivity of crystalline enstatite. 

3.1.2. NGC6537 

The results for NGC6537 are shown in Fig. ||. The temper- 
atures found for the forsterite (75 K) and enstatite (60 K) 
in NGC6537 are among the lowest found in our sample. 
Note that if an extra dust component significantly con- 
tributes to the 40 micron complex, the temperature of 
enstatite will be higher (and its mass lower) than what 
has been determined here. 

The spectral energy distribution of the complete spec- 
trum is too broad to be fitted by a single temperature dust 
component. 

3.1.3. NGC6302 

The continuum subtracted spectrum of NGC 6302 and its 
good fit are shown in Fig. ^. 

Molster et al. (2001b) used the same method as used 
in this paper, and therefore found the same temperatures. 
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Fig. 7. A fit (dotted line) to the continuum subtracted spec- 
trum (solid line) of NGC6302. T/ = 65 K and Te = 70 K. 



As for NGC6537, it was not possible to fit the spectral 
energy distribution with a single temperature dust com- 
ponent. Molster et al. (2001b) attribute the broad energy 
distribution to the presence of a population of large grains, 
which mainly contribute to the long wavelength side. The 
presence of this population of large grains is indicated by 
the gentle slope of the spectrum up to mm wavelengths 
(Hoare et al. 1992). 

The temperature found for the enstatite and forsterite, 
respectively 65 and 70 K, are similar to the temperature 
of NGC6537, which in many aspects looks very similar to 
NGC6302. Kemper et al. (2001) assumed two temperature 
regimes: a cold one from 30 to 60 K, and a warm one 
from 100 to 118 K. Both components contain forsterite and 
enstatite. Our results, giving a temperature somewhere in 
between those two regimes, is in agreement with theirs, 
although the exact comparison is somewhat difficult. 

3.1.4. MWC922 

The fit to the continuum subtracted spectrum of MWC922 
is one of the best we have (see Fig ||). Especially the 40 
micron complex is very well reproduced by our model, 
indicating that the 50% clino- and 50% ortho-enstatite 
are the right proportions for this object. At A < 16/zm 
the spectrum of MWC922 is dominated by PAH-features 
which were not incorporated in the fitting procedure. 

3.1.5. AC Her 

A model with cool dust fits the long wavelength part 
(> 2G^m) of the AC Her spectrum (dotted hne in Fig. 
However, the short wavelength features indicate the pres- 
ence of a dust component with a much higher temperature. 
The temperature of this material is not well constrained. 
In Fig. ^ we show a fit of 700 K (dashed line in Fig. 
but a similar fit could be derived with a temperature sev- 
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Fig. 8. A fit (dotted line) to the continuum subtracted spec- 
trum (solid line) of MWC922. T/ = 90 K and Te = 100 K. 
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Fig. 9. A fit (dotted line) to the continuum subtracted spec- 
trum (solid line) of AC Her. Tf = 100 K and Te = 90 K 
(dotted line). The dashed line is a 700 K (for both forsterite 
and enstatite) model fit. 

eral hundreds degrees Kelvin higher or lower. Therefore it 
is impossible to give a reliable mass estimate for this hot 
component. 

In our modelling we only assumed a single tempera- 
ture. Based on the necessity of (at least) two different tem- 
peratures, the existence of a temperature gradient seems 
more likely. It is interesting to note that the overall spec- 
trum of AC Her is very similar to that of comet Hale Bopp 
(Molster et al. 1999a) where we know that the dust is lo- 
cated in one place. Temperature differences found in the 
grains around this comet must therefore originate from 
the grain size differences. Small grains can account for 
the high temperature dust emission, while bigger grains 
are responsible for the low temperature dust emission. 
Such a scenario might also be possible for AC Her, which 
would imply that the dust might not have to be so close to 
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Fig. 10. A fit (dotted line) to the continuum and amor- 
phous silicate subtracted spectrum (solid line) of HD45677. 
Tf = 140 K and Te = 140 K. 



the star as previously thought (e.g. Alcolea & Bujarrabal 
1991). Jura et al. (2000) found a disk Hke structure for this 
object, which supports the above mentioned scenario. A 
full radiative transfer model fit would be necessary to com- 
pletely understand the dust distribution around AC Her, 
but that is beyond the scope of this paper. 

3.1.6. HD45677 

From the continuum subtracted spectrum of HD45677 
we first removed the broad amorphous silicate features 
(Fig. |l^). We cannot exclude that we also removed part 
of the crystalline silicate features in the 18 micron com- 
plex in this way. This does not influence our results since 
these are mainly based on the 23, 28, 33 and 40 micron 
complexes. To fit the spectrum of HD45677 we ignored 
the strength of the 19.5 micron feature, which is severely 
overestimated in our resulting fit. If we would have fit- 
ted the 19.5 and 40 micron features simultaneously, the 
28 micron complex would have been severely underesti- 
mated. Likewise, attempts to fit the 18 and 28 micron 
complex together will result in a severely overestimated 40 
micron complex, and also the fits to the 23 and 33 micron 
complexes will become worse. It is unlikely that this dis- 
crepancy can be fully explained by the subtraction of the 
amorphous silicates. Because this is not the only source 
with this problem, we leave this for future research. 

Malfait (1999) also studied this star. He modelled this 
object with a radiative transfer code. HD45677 could only 
be modelled with a 2 component dust shell, consisting of 
a hot shell, responsible for the main part of the flux up to 
20 /xm and a cool component which is the main contributor 
to the crystalline silicates features. Due to the method we 
use here, our temperature estimate is based on this cool 
component. Malfait finds a temperature between 250 and 
50 K for this cool component. Unfortunately this is not 
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Fig. 11. A fit (dotted line) to the continuum subtracted spec- 
trum (solid line) of 89 Her. Tf = 110 K and = 100 K. 



specified for the different components separately, so we 
can only say that our temperature estimates do agree with 
this temperature range. 

The predicted strength of the crystalline silicate fea- 
tures in the 10 micron complex is underestimated. Since 
the strength of the amorphous silicate band at 10 /Ltm is 
uncertain, errors in the estimate of its contribution affect 
the strength of the crystalline silicate bands at these wave- 
lengths and we did not attempt to fit the hot crystalline 
silicate compounds. 

3.1.7. 89 Her 

Before we fitted the continuum subtracted spectrum of 
89 Her, we first subtracted a broad feature below the 26 
to 45 /um region (Fig. [Tl] ). This feature is also seen in 
HD44179 and probably AFGL 4106 and discussed in Pa- 
per II. 

The continuum subtracted spectrum of 89 Her is quite 
noisy at the longer wavelengths, which makes the fits not 
as well constrained as in other stars. Also in this star 
warmer grains are necessary to explain the crystalline sili- 
cate structure found on top of the amorphous silicate fea- 
ture in the 10 micron complex. Again, problems in the 
separation of the crystalline and amorphous silicates kept 
us from fitting this feature. Based on the CO observations 
and the near-IR excess, it was argued in Paper I, that 
there must be dust with different temperatures, likely a 
temperature gradient, around 89 Her. 

3.1.8. MWC300 

Although we have argued in this paper that the strength 
of the 19.5 micron feature is difficult to model correctly, 
we decided, because of the problems in the 28 micron com- 
plex to constrain the enstatite by the 19.5 micron feature 
in MWC300 (see Fig. |l|). If we would have fitted the 
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Fig. 12. A fit (dotted line) to the continuum subtracted spec- 
trum (solid line) of MWC300. T/ = 90 K and T, = 90 K. 




Fig. 13. A fit (dotted line) to the continuum subtracted spec- 
trum (solid line) of HD44179. T/ = 135 K and T, = 135 K. 
Below 15 /xm the spectrum is dominated by PAH features. 



Fig. 14. A fit (dotted line) to the continuum subtracted spec- 
trum (solid line) of HD161796. Tf = 105 K and Te = 80 K. We 
also included crystalline water ice with a temperature of 70 K 



3.1.9. HD44179 

As for 89 Her, we removed the very broad feature in the 
26 to 45 /im range in the continuum subtracted spectrum 
of HD44179. The 18 micron complex seems to contain a 
contribution from amorphous sihcates which was also re- 
moved. The result can be found in Fig. |l^. It should be 
noted that a change in the subtraction of the broad 18 
micron amorphous silicate feature, whose properties are 
not well determined, may change the derived temperature 
by more than the typical fitting error of 10 K. 

The derived continuum temperature of 120 K is some- 
what uncertain due to the complex nature of the source 
(Waters et al. 1998). Since the short wavelength part of the 
continuum is formed by C-rich grains, we based our fits on 
the long wavelength range. This will likely underestimate 
the temperature of the continuum. Taken all these un- 
certainties into account, the crystalline silicates may have 
an equal or even lower temperature than the amorphous 
silicates, as is found in the other stars. 



strength of the 28 micron complex, and ignored the 19.5 
micron feature, we would have derived a temperature of 
roughly 150 K, which is much larger than the forsterite 
temperature. It would also predict prominent features at 
the shorter wavelengths, which were not seen. The use of 
the 19.5 micron feature to constrain enstatite resulted in 
a similar temperature for the forsterite and enstatite dust 
species. This result, together with the reasonable fit at the 
40 micron complex makes us confident in our approach for 
this star. We note however, that because of the problems 
with the 28 micron feature the values for enstatite are 
poorly constrained. The source of the extra flux in the 28 
micron feature is unknown. 



3.1.10. HD161796 

The 40 micron complex in HD 161 796 is dominated by 
crystalline water ice. In order to fit this spectrum (Fig. Q 
we therefore added crystalline water ice to the spectrum 
(Smith et al. 1994). In general most features are reason- 
ably well reproduced. The poor fit around 40 /xm suggests 
that an underlying weak broad component contributes. 
Possibly, this is amorphous water ice. 

Within the errors, the crystalline and amorphous tem- 
peratures are the same. The peak wavelength of the crys- 
talline water ice feature in our model is slightly offset 
from our ISO spectrum. This may be a temperature ef- 
fect (Smith et al. 1994), reflecting the sensitivity of the 
emissivities to temperature. 
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Fig. 15. A fit (dotted line) to the continuum subtracted spec- 
trum (solid line) of HD179821. T/ = 75 K and Te = 65 K. 
To fit the 40 and 60 micron complex, we included crystalline 
water ice with a temperature of 45 K. Note, that we made no 
attempt to fit the amorphous silicate bands in the continuum 
subtracted spectrum 
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Fig. 16. A fit (dotted line) to the continuum subtracted spec- 
trum (solid line) of AFGL4106. Tf = 100 K and = 80 K. 
Note, that we made no attempt to fit the amorphous silicate 
bands in the continuum subtracted spectrum 



3.1.11. HD179821 

HD179821 shows prominent crystalline water ice features 
in the 40 and 60 micron complex. Therefore, we also added 
crystalline water ice in this model fit (Fig. |l^). The crys- 
talline silicate features are rather cold, 65 K for the en- 
statite and 75 K for the forsterite, therefore no detectable 
features are expected in the 10 micron region. 

3.1.12. AFGL4106 

For the fit to the continuum subtracted spectrum of 
AFGL4106 (Fig. ^6|) we changed the ratio of clino versus 



ortho enstatite to 1 : 3. Otherwise the 44.7 micron feature 
would have appeared too strong in the model spectrum. 
This feature is located in the wing of the 43.0 micron fea- 
ture. The 43.0 micron feature is not well reproduced by 
our simple model, resulting in an offset of the 44.7 micron 
feature. 

A broad feature between 30 and 45 fj,m seems present. 
It is not known whether this is the same feature as in 89 
Her and HD44179. It seems to peak at a longer wavelength 
(38 /im) than in the other two sources. Subtracting this 
feature, would increase both the temperature of forsterite 
and enstatite. The 40.5 /im feature is very strong, even 
if one removes this broad feature, indicating that other 
dust species are present. No attempt has been made to 
fit the substructure in the 18 micron band in view of the 
large and uncertain amorphous silicate contribution. In 
this star no detectable 10 micron structure is expected 
from the crystalline silicates, which is in agreement with 
the observations. 

A full radiative transfer calculation was made for this 
source by Molster et al. (1999b). The enstatite abundance 
derived by them was based on the (wrong) assumption 
that the 32.8 micron feature was due to ortho-enstatite, 
and is therefore difficult to compare. This is related to 
their use of another datasets for forsterite, which has a 
different intrinsic ratio of the 23.7 and 33.6 micron fea- 
tures. The result is a lower temperature and therefore 
higher mass for forsterite in their calculations. Also, for 
the amorphous olivines they used a different dataset, how- 
ever in this case it resulted in a higher temperature and 
therefore lower mass. The combination of those two ef- 
fects explains the low amorphous silicate to forsterite ratio 
found by them compared to this study. 



3.1.13. NML Cyg 

We could not reliably fit enstatite to the spectrum of 
NML Cyg since a part of the 28 micron complex was miss- 
ing, due to problems with the data reduction (see Paper I). 
We were able to fit forsterite to the continuum subtracted 
spectrum and found a temperature of about 150 K (not 
shown) . 

3.1.14. IRC-K10420 

Also for IRC+ 10420 we only fitted the forsterite compo- 
nent which appeared to be about 90 K (not shown). En- 
statite could not be fitted due to the lack of data in the 
28 micron region. We checked the 40 micron complex and, 
if enstatite would have been fitted, an ortho-enstatite to 
clino-enstatite ratio of 2 : 1 would probably give the best 
fit to the strength of the 44.7 micron feature. 
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Fig. 17. The temperature of the forsterite grains versus the 
wavelength position of the 69 micron band. The circles are 
outflow sources and the diamonds are disk sources. 



4. Correlations 

From the modelling, we were able to derive values for the 
temperature and the relative contributions to the total 
amount of dust mass by forsterite and enstatite. Corre- 
lating these values can give interesting insight in the dust 
inventory as we will show in this section. Since in Paper 
I and II it was found that there is a significant difference 
between the spectra of disk and of outflow sources, we 
will separate them. In all plots in this section the outflow 
sources are represented by a circle while the disk sources 
are represented by diamonds. 



Temperature dependence of laboratory spectra: peak 
width and position 



In Sec |2.1| we showed a relation between the peak posi- 
tion and FWHM of the 69 micron forsterite band and the 
temperature of the forsterite grains. Since we have now 
determined the temperature of the forsterite grains, we 
checked this relation in our data (Fig. |l^. There is no 
clear correlation between these two quantities. A possible 
explanation for this scatter behaviour might be the fact 
that the temperature of the forsterite is determined from 
bands in the 20 to 40 /im range. The strength and width of 
the 69 micron feature may be dominated by much cooler 
dust. The predicted strength from our simple model fits, 
which is lower than the strength in our ISO spectra (see 
e.g. HD179821), supports this statement. 

4-. 2. Other temperature trends 



In Figure 18 we compare the temperature of the enstatite 
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Fig. 18. The temperature of the forsterite grains versus the 
temperature of the enstatite grains. The diamonds are the disk 
and the circles are the outflow sources. The dashed line repre- 
sents equal temperatures for the forsterite and enstatite grains. 
The errors are typically in the order of 10 K 
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and forsterite, derived from our simple model fits. For the 
disk sources the enstatite and forsterite grains seem to 
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Fig. 19. The temperature of the forsterite grains versus the 
temperature of the amorphous silicate grains. The dashed line 
represents equal temperatures for the forsterite and amorphous 
silicate grains. The diamonds correspond to the disk sources, 
while the circles correspond to the outflow sources. 



have an equal temperature, while in the 3 outflow sources 
the forsterite seems warmer. 

In Fig. |l| we compare the forsterite temperature and 
the amorphous silicate temperature. We first note that, in 
general, the crystalline forsterite grains are colder than the 
underlying continuum consisting of amorphous silicates. 
The difference in temperature between the amorphous and 
crystalline silicates has probably to do with the difference 
in chemical structure. The crystalline silicates are Mg-rich 
(see this paper and Paper II), while the amorphous sili- 
cates must contain metals to explain their higher near- 
infrared absorptivity. Possibly the reaction of the Mg-rich 
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Fig. 20. The mass of the forsterite grains versus the mass 
of the enstatite grains. We have assumed that they have the 
same size and shape distribution. The dashed line represents 
equal masses for the forsterite and enstatite grains. The outflow 
sources are marked as a circle while the disk sources are marked 
with a diamond. 

crystalline silicates with gaseous iron may proceed in these 
outflows at temperatures well below the glass temperature 
leading simultaneously to amorphous and dirty silicates 
(Tielens et al. 1998). 

For some stars like IRC+10420, MWC922 and 
IRAS09425-6040 the difference in temperature excludes 
the possibility of one grain population, which is partially 
crystalline. The crystalline and amorphous silicates must 
form two separate grain populations, which are not in 
thermal contact. From this simple model we cannot say 
that they are also spatially separated, but radiative trans- 
fer modelling (Molster et al. 1999b; 2001a) indicates that 
this is not necessary. Finally, we note that there is no obvi- 
ous separation between the disk and the outflow sources. 
Because of the analogy between the different sources, we 
expect that this grain segregation is valid for all sources 
in our sample. 

4-. 3. Abundance trends 

From our modelling we were also able to derive an abun- 
dance ratio for enstatite and forsterite (Fig. ^). In all 
sources, except NGC6302, the enstatite is more abundant 
than the forsterite, on average a factor 3-4. Again there 
seems to be a difference between the disk and outflow 
sources. The outflow sources have, on average, a higher 
enstatite over forsterite ratio than the disk sources. We 
note however, that the number of outflow sources is low 
and more data is required to confirm this trend. 

If crystallization of the amorphous grains is a ther- 
mally driven process, it is expected that warmer dust 
shells would show a higher crystalline silicate fraction. 
This effect would be most clear in the disks sources, where 



the dust stays at roughly the same place. In the out- 
flow sources, the temperature of the optically thin dust 
shell is more an indication of the age of the dust. Cooler 
dust would be further out and therefore older, in which 
case an increase of the crystallinity with a lowering of 
the temperature may even be possible. We investigated 
the crystalline over amorphous silicates abundance ratio 
with the temperature of the amorphous silicates. We ex- 
cluded AC Her, 89 Her and HD45677 from this analysis 
because the continuum and feature temperatures derived 
likely do not reflect particles with the same size distri- 
bution and/or location around these objects, resulting in 
unrealistic mass ratios. No clear relation could be found 
between the temperature of the amorphous silicates ver- 
sus the crystallinity, neither for the disk nor for the out- 
flow sources. The crystallization process is either a relic 
of the past or it is not thermally driven. This is not too 
surprising. Thermally driven crystallization requires tem- 
peratures of « lOOOK (Hallenbeck & Nuth, 1998), which 
is well above the temperature of either the crystalline or 
amorphous silicates in our sources. 

4-4- Crystallinity versus stellar flux ratio 

It has been suggested by Waters et al. (1999) that the 
enstatite over forsterite ratio increases with luminosity of 
the star. They based their conclusions on the ratio of the 
32.97 and 33.6 micron bands. They attributed the former 
band to enstatite, and the latter one to forsterite. Unfor- 
tunately, the 32.97 micron band is strongly affected by 
instrumental effects (see Paper II) , which are most promi- 
nent in the brightest sources. These are usually also the 
intrinsic most luminous sources, such as IRC4-10420 and 
NML Cyg. Since it is likely that in high flux sources grain 
formation will result in different ratios of the dust species 
we still decided to investigated this scenario. We could find 
no trend for the stars in our sample, which had a reliable 
distance (and thus flux) estimate. 

5. Discussion 

Most features can be explained with forsterite and en- 
statite. However laboratory spectra of these species often 
produce too broad bands. Most laboratory spectra were 
measured at room temperature, while the dust species 
around our stars are in the order of 80 K. Lowering the 
temperature can narrow the dust features, and this might 
be the natural solution for the width problem. It can 
also solve another problem. The 69.0 micron feature has 
been attributed to forsterite. It is very sensitive to the 
Fe/(Mg-|-Fe) ratio in olivine, and one of the key features 
for the determination of the Fe content in crystalline sil- 
icates. However in most lab spectra it peaks at 69.7 /Ltm. 
Recent laboratory measurements show that this feature 
shifts to shorter wavelengths when the sample is cooled 
and will be at w 69.1 /im at 77 K and even at « 68.8^m 
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when measured at 4 K (Bowey et al. 2000). At the same 
time, the feature narrows (cf. Fig ^ . So it seems that the 
temperature dependence is able to solve also this problem. 

From the correlation study in this chapter a clear dif- 
ference between the disk and the outflow sources emerges. 
The differences found in Paper I and Paper II were based 
on the shape and strength of the different solid state 
features in the spectrum. These differences can now be 
traced back to differences in temperature and chemical 
composition of the circumstellar dust. The outflow sources 
seem more abundant in cold enstatite grains than the disk 
sources, although more observations would help to quan- 
tify this better. There are also indications that in the high 
luminosity outflow sources the ortho- over clino-enstatite 
abundance ratio is larger than unity (see AFGL4106 and 
IRC-l-10420), while in the other (disk) sources nice fits 
were obtained with a ratio of unity. Because of the large 
crystalline water-ice component in the 40 micron complex 
of IID161796, it is difficult to quantify its ortho- over clino- 
enstatite ratio. Therefore we cannot exclude that this dif- 
ference we see is related to the more massive nature of the 
central stars of AFGL4106 and IRC+10420 rather than 
the disk-outflow character. 

The difference between the disk and the outflow 
sources might be related to differences in conditions dur- 
ing the condensation of the grains out of the hot gas, 
and/or to differences in the conditions since the forma- 
tion of the dust particles. Smyth (1974) found that clino- 
enstatite slowely inverts to ortho-enstatite between 920 
and 1220 K. Dust particles will stay longer within this 
temperature range in the case of massive stars than for 
low mass stars. This might explain the overabundance 
of ortho-enstatite in IRC-H10420 and AFGL4106 with re- 
spect to the other stars. Another possible difference during 
the formation of these grains could be the amount of ra- 
diation pressure exerted on the particles. It is expected to 
be lower in the disk sources, otherwise the disk is likely to 
be blown away. Also the presence of a companion -most if 
not all disk sources are a binary system- might influence 
the dust forming process. If, on the other hand, we as- 
sume that the initial (forming) conditions are similar for 
the disk and outflow sources, than the differences in the 
conditions after formation must be dominant. Time is an 
obvious differences. Also a relation with grain coagulation 
-the disks contain large grains- cannot completely be ruled 
out. 

We note that clino-enstatite is often found in me- 
teorites on earth. The process responsible for the over- 
abundance of clino-enstatite in meteorites in the disks 
of young stars might be the same as the one in the 
disks around evolved stars. We also want to note that in 
chondrites or primitive meteorites, i.e the less processed 
ones, clino-pyroxenes are less abundant than the ortho- 
pyroxenes. This similarity between the disks around young 
and around evolved stars, is new evidence to the hypoth- 
esis that the circumstances and processes in disks around 



young stars are very similar to these in disks around 
evolved stars, although their origin is quite different. 

The equal temperature of enstatite and forsterite in 
the disk sources, is in principle compatible with the as- 
sumption that enstatite and forsterite are present as a 
composite grain in the disk sources. However, it is equally 
well possible that their individual equilibrium tempera- 
tures are similar, in which case they can still be individ- 
ual grains. New uninterupted laboratory measurements of 
the optical properties of both materials, especially in the 
wavelength range were they absorb the stellar light and 
emit the thermal radiation, would allow us to make these 
kind of calculations and settle this point. The different 
temperatures found for forsterite and enstatite in the out- 
flow sources suggests two separate grain populations in 
these environments. 

6. Conclusions 

We can summarize the main results of this study as fol- 
lows: 

1 The ISO spectra could be reasonably well fitted with 
laboratory spectra of forsterite and enstatite. 

2 The models underestimate the flux at 18, 29.6, 30.6, 48 
^m and sometimes at 40.5 /xm, which is an indication 
for the presence of (an)other dust component(s). No 
convincing identiflcation could be made yet. Diopside 
does have features at most of these wavelengths, but 
also strong features at others which are weak or absent 
in the ISO-spectra. 

3 The 19.5 micron feature is often overestimated by our 
model spectra. No explanation is yet known for this 
phenomenon, but it should be noted that in the full ra- 
diative transfer modelling it appeared to be much less 
of a problem. This might indicate that optical depth ef- 
fects play a role. Also the calculation of the absorption 
coefficients from the constants instead of the absorp- 
tion coefficients from laboratory particles might lead 
to differences. 

4 The band width of the laboratory data is larger than in 
our ISO spectra. This difference is likely a temperature 
effect, and might be used as an independent tempera- 
ture indicator. Especially the 69.0 micron band is very 
suitable for this analysis. 

5 The temperature of the forsterite and enstatite grains 
are (almost) similar for the disk sources, while the 
forsterite is slightly warmer in the outflow sources. 
This would imply that the forsterite and/or enstatite 
grains differ slightly in the disk and outflow sources. It 
is not clear whether this difference is due to a differ- 
ent formation process, or due to a different dust pro- 
cess history after the grain formation. Since this trend 
is only based on 3 sources more data is required to 
confirm the difference between the dust and outflow 
sources. 
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6 The crystalline silicates are colder than the amorphous 
silicates. This is probably due to the difference in 
chemical composition. No Fe is present in the crys- 
talline silicates, while in the amorphous silicates it is 
expected to explain the higher absorptivity. This dif- 
ference in temperature also implies that the crystalline 
and amorphous grains are two distinct grain popula- 
tions. 

7 Enstatite is on average a factor 3-4 more abundant 
than forsteritc in our sources. There are indications 
that the enstatite over forsterite ratio in the outflow 
sources is higher than in the disk sources. 

8 In the low luminosity sources the spectra were well fit- 
ted with an equal amount of ortho- and clino-enstatite, 
while in the two high luminosity sources more ortho- 
enstatite seems to be present. 

9 No correlation could be found between the crystallinity 
and the temperature of the dust. Also the luminosity 
of the stars does not seem to be correlated with the 
enstatite over forsterite ratio. 

10 These simple model fits already give a good insight 
in the dust around our stars. In Paper I the shape of 
the features naturally separated the disk and outfiow 
sources. In this study the differences between these two 
categories became again evident. 
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